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CODSTRAN: COMPOSITE DURABILITY STRUCTURAL ANALYSIS 


C. C. Chamis and G T. Smith 
NASA-Lewls Research Center 


ABSTRACT 

CODSTRAN (Composite Durability STRuctural ANalysis) la an Integrated computer program 
being developed for the prediction of defect growth and fracture of composite structures 
subjected to service loads and environments. CODSTRAN is briefly described herein with 
respect to organization, capabilities and present status. Application of CODSTRAN current 
capability to a flat composite laminate with a center slit which was subjected to axial 
tension loading predicted defect growth which is in good agreement with C-scan ultrasonic 
test records. 


INTRODUCTION 

Assessing and establishing the structural integrity and durability of composite 
structural components in service environments requires analysis capabilities which can re- 
liably predict defect growth and fracture in these components. The defect in these com- 
ponents may be ar initial flaw or it may be caused by stress concentrations. For quanti- 
tative predictions of defect growth and fracture in composite structural components the 
analysis capabilities should account for: complex component shapes, complex loading condi- 

tions (such as static, cyclic, dynamic and environmental), and material and geometric non- 
linearities. In addition, the analysis capabilities should be relatively easy to use and 
update. These types of capabilltir , by necessity, will be in the form of Integrated 
computer programs. An integrated .iputer capability with many of the above features is 
under development at NASA Lewis Research Center (LeRC) under the acronym CODSTRAN (Com- 
posite Durability STRuctural ANalysis). The objective of this paper is to provide a brief 
description of CODSTRAN with respect to its organization, capabilities, present status, 
and applications to date. 


CODSTRAN 

CODSTRAN is an integrated computer program designed to predict defect growth and 
fracture of composite structural components subjected to service load and environmental 
conditions. CODSTRAN is modular and open ended. Present modules may be replaced as im- 
proved ones becomes available and new ones can be added as they are developed. The organ- 
ization (modular structure) of CODSTRAN, the function of its various modules and its capa- 
bility are described in some detail in this section. 


\ 
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ORGANIZATION 


The organization of CODS TRAN la shown schematically In the flow chart, figure 1. As 
can be seen In this figure CODS TRAN consists of seven major modules: Executive, Input, 

Output, Analysis, Composite Mechanics, Fracture Crlrerla, and Life Prediction. Communica- 
tion between modules Is through the Executive module. A brief description of the modules 

fol lows . 

Executive module . The Executive module contains all the execution Instructions of 
CODSTRAN and the communication links between the different modules. The user controls the 
program through the executive module. The logic flow of the execution Instructions for 
defect growth and fracture, using Input (I/p), Output (0/p), the composite mechanics In 
MFCA (Multilayer Filamentary Composite Analysis Code (ref. 1)) and the structural/streas 
analysis In NASTRAN (ref. 2), Is shown In figure 2. Note that the blocks are numbered In 
the upper left corner for convenience of reference. Each block In this diagram covers one 
or several execution Instructions. For example: laminate fracture is Identified In block 9 

via the stability check; ply failure In block 20; failure mode In bloc u . 21; lnterply de- 
lamination In block 23; and defect growth In block 29. The crack opening displacement Is 
calculated In block 12. The grid (node) coordinates and element stresses are updated In 
block 16. The load redistribution Is accomplished by bypassing block 30 when defect growth 
has occurred . 

Input module The Input module (I/p) controls all the input Including component geom- 
etry, defect geometry, composite type, laminate configuration, material properties, fracture 
data, fatigue data, and service environment. The input Is expedited by using special pur- 
pose preprocessors (such as mesh generators) and computer resident data banks where compos- 
ite data are compiled as they become available. The I/p module has provisions to activate 
available preprocessors and computer resident data banks. For example, the computer resi- 
dent data bank for use with the MFCA Is also available to CODSTRAN through the Input module. 
This data bank currently has properties for 14 different fibers, 10 different matrices, 

10 different unidirectional composite systems, and 3 intraply hybrids. 

Output module . The Output module (0/p) has the operational Instructions for printing 
out the analysis results such as crack opening displacement, failed plies, fracture mode 
inducing failure, and defect growth. The output data are printed out In convenient tabular 
or graphical form for ease of Interpretation. Additional and/or Intermediate data may be 
printed out at the user's request. All the 0/p features of NASTRAN are available. 

Analysis modules . The Analysis modules planned for CODSTRAN consist of the following: 

(1) approximate analysis modules based mainly on the mechanics of materials methods; 

(2) fracture mechanics modules based on the 1/^r (r is the distance from the crack tip) 
near-fleld stress variation for slits and comparable ones for holes; (3) semi-emplrical 
equations obtained by fitting strain gage data for the near field strain variation; (4) 

(4) special finite elements (hybrid singular element, ref. 3); and (5) NASTRAN (ref. 2). 
NASTRAN can be used for both the near- and far-fleld solution of large structures or 
structural components while the methods listed. Items (2) to (4), are used only locally. 

Composite mechanics module . The Composite Mechanics module consists mainly of the 
MFCA Computer Code (ref. 1). This code has composite mlcromechanlcs , composite macro- 
mechanics, combined stress failure criteria, linear and nonlinear laminate theories, and 
lnterply (interlaminar) failure (delamination) criteria. In addition, provisions are made 
for: (1) residual stress computations (ref. 4); (2) environmental (temperature and mois- 
ture) effect (ref. 5); (3) In situ ply strengths (ref. 6); and (4) free edge delamlnatlon 
due to interlaminar stresses using approximate relationships (ref. 7). The Composite 
Mechanics module Is called several times during the analysis process. It is called to gen- 
erate laminate properties for finite element analysis, compute ply and lnterply stresses, 
compute residual stresses, and check for ply and lnterply failures. Referring to figure 2, 
execution Instructions calling for composite mechanics Include blocks 3, 4, 18 to 24, 27 
and 28. 
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Fracture Criteria module . The Fracture Criteria module la designed to accommodate 
atatlc cyclic and dynamic fracture The criteria for static fracture are of two forms: 
one for ply level fracture and the other for laminate fracture. The ply level fracture 
criteria are the two available In MFCA (ref. 1). One la based on a modified distortion en- 
ergy principle and the other on a general quadratic surface fit. The interply fracture 
criteria used Is that which Is available in MFCA and Is based on limiting the relative ro- 
tation of adjacent plies. The laminate level fracture criteria consist of the Inherent 
flaw (refs. 9 and 10). In addition, provisions are available for several other fracture 
criteria such as the minimum strain energy density and combinations of linear elastic frac- 
ture mechanics modes I, IT, and III (ref. 11). 

Cyclic and dynamic fracture criteria are not available as yet. Hlgh-strain-rate ex- 
perimental data currently being generated under a NASA contract will be used to develop ap- 
proximate dynamic fracture criteria. Data from the literature are currently being examined 
to Identify relevant cyclic failure criteria. 

Life-Prediction module . The Life Prediction module will eventually provide the formal- 
isms to predict the life and/or durability of a composite part under sustained load, fatigue, 
and service environments. No life prediction formalisms of this type are available as yet. 
The Wear-Out module (ref. 12) appears to be a step In this direction. Data currently gen- 
erated under another NASA LeRC program will be used to develop life prediction criteria. 


Capabilities 

The functional features of the various modules will eventually provide CODSTRAN with 
the following capabilities: 

1. Durability assessment of large structures and complex structural parts from com- 
posites , 

2. Durability of components with regular and irregular defect geometry, 

3. Accurate predictions of stress states near defects and discontinuities, 

U . Structural response due to static, cyclic, transient impact and thermal loads, 

3. Evaluation of defects in all types of fibrous composites and hybrids, 

6. Assessment of environmental (temperature and moisture) effects on fibrous compos- 
ites with defects, 

7. Assessment of geometry and material nonlinearities on defect growth and fracture, 

8. Criteria for static, cyclic and dynamic fracture, 

9. Ply, Interply and free-edge interlaminar failure criteria, 

10. Predicting laminate level fracture when defects are present, 

11. Predicting residual stress effects on defect growth, 

12. Evaluating in situ ply strength effects on defect growth. 

It is belived at this time that these capabilities are essential for reliably pre- 
dicting the structural Integrity of fibrous composite components and for assessing the as- 
sociated durability. 
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Present Status 


CODSTRAN is being developed at LeRC on the UNIVAC 1110 computer. The Executive module 
is almost completed for communication between the I/p, 0/p, NASTRAN and MFCA modules Some 
debugging, storage file and system associated difficulties exist in updating the solution 
from one load increment to the next after several cycles have been completed. At this 
stage of development, CODSTRAN lias the following analysis capabilities: 

1. Complex structural components, flaw geometry and load conditions - handled via 
NASTRAN available capabilities, 

2. Composite mechanics - available in MFCA, 

3. Static failures - ply-by ply and tntraply through MFCA, 

4. All types of fibrous composites and several types of hybrids - handled through MFCA, 

5. Lamination residual stresses - computed by MFCA, 

6. Geometry and material nonllncarlties - incorporate at the load Increment following 
defect growth through NASTRAN and MFCA. 

Even at this early stage of the development CODSTRAN has considers! .e capability for 
the prediction of defect growth and fracture in fibrous composite structures. 


APPLICATION 

CODSTRAN, with its present capability, has been applied to predict the defect growth 
in a flat specimen with a center through-s 1 it and subjected to tension. The data used are 
from reference (13). The specimen geometry, composite system, and laminate configuration 
are shown in figure 3. This specimen has a width to crack-length ratio of about 5. The 
finite element model Is shown In figure 4. The whole specimen was modeled in order to 
minimize the effects of displacement constraints that would have been required had advan- 
tage of symmetry been taken. The model consists of 530 plate elements (NASTRAN Library 
CQL'AD2 and TRIAG2), 442 nodes (grid points) and 870 degrees of freedom (DOF) Two degrees 

of freedom per node are allowed for lnplane loadings. 

The variation of the stress o predicted by CODSTRAN is plotted as a function of the 
distance "r" from the crack tip in figure 5. Two significant stress field characteristics 
in figure 5 are: (1) The normal stress varies inversely as the square root of the distance 

( \lxfr ) from up to r/a - 0.7 (2a is the crack length which corresponds to a distance of 

about two laminate thicknesses (0.4a) from the crack tip; and (2) the normal stress Oy 
approaches thr gross stress at a distance of about four laminate thicknesses from the 

crack tip. Both of these observations illustrate thepolnt that the stress concentration 
predominates in a relatively small region (about one laminate thickness) near the defect 
boundary . 

CODSTRAN results of partial defect growth at 1000 pound intervals are shown in fig- 
ures 6(a) to (h). In these figures, the initial crack is shown by the darkened squares. 
Element ply failures (locally damaged areas) are shown by the element outline, (open 
elements). In these elements any one or combinations of the following ply failures have 
occurred: transverse tension, intralaminar shear, fiber fractures, and interply delamina- 

tion Furthermore, these failures may have occurred in one or more plies or interply 
layers, but not in all the plies. A ply within the element has failed completely when it 
cannot carry additional load in any of the three fracture modes: longitudinal, transverse 
and intralaminar shear. An clement has failed completely when all of its plies have 
failed completely. Complete element failure indicates defect growth (crack opening) and 
is depicted with the darkened elements in the computer plots. 
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Aa can be aeen In figure 6(b) the two thousand pound load haa cauaed ply fallurea In 
eight elementa at the crack enda . Flgurea 6(c) to (f) ahow how ply fallurea progreaa away 
from the crack. The nonaymmetrlea aeen In theae flgurea are due to the tolerancea (* 5 per- 
cent) allowed In the failure criteria. Figure 6(g) ahowa four complete element fallurea 
at eight thousand pounds. Note the extensive local ply fallurea shown In figure 6(h). The 
results In these figures taken collectively show that CODSTRAK and, therefore, the Inte- 
grated theory It represents, predicts progressive failure and defect growth In composite 
laminates with defects. 

C-Scan ultrasonic test records of the laminate prior to load (ref. 13) and at a load 
of about one-half the fracture load (11 000 lb) are shown In figure 7. This corresponds 
approximately to the damage growth predicted by CODSTRAN shown In figure 6(e). Figures 7(b) 
and 6(e) are shown In figures 8(a) and (b), respectively for comparison purposes. The 
undamaged part In figure 8(b) has been darkened to correspond to the C-Scan record 
(fig. 8(a)). It can be seen that the predicted results are In good qualitative agreement 
and are In reasonable quantitative agreement with the C-scan records. The predicted damage 
results (fig. 8(b)) are narrower and longer compared to the C-scan records (fig. 8(a)). 

The physical Interpretation of this Is that the predictions Indicate less delamlnatlon 
damage and more transverse failure compared to the C-scan. The C-scan would not pick up 
the transverse failure readily since the record was made when the specimen was not under 
load. 

The crack opening displacements predicted by CODSTRAN are also compared with those 
that were measured as shown In figure 9. As can be seen In this figure, CODSTRAN predicts 
somewhat smaller crack opening displacements compared to measured data. Though this may be 
considered as a reasonable agreement for the Initial evaluation of CODSTRAN, It Is clear 
that C0DS T RAN needs additional refinements and tuning for Improved quantitative correlations. 

The results presented and discussed In this section demonstrate, In a limited way, 
some of the computational features, the potential, and the development status of CODSTRAN. 


SUMMARY OF RESULTS 

This paper presents a brief description of CODSTRAN (Composite Durability STRuctural 
ANalysis), an Integrated computer program under development at NASA Lewis Research Center. 

It is designed to predict defect growth and fracture of composite structures subjected to 
service loads and environments. When completed, CODSTRAN will account for geometry and 
material nonlinearities, environmental effects as well as static, cyclic and dynamic frac- 
ture. At thepresent state of development, CODSTRAN-predicted defect growth Is in reasonable 
agreement with that observed in flat composite laminate:* with center silts and subjected to 
uniaxial tension. 
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Figure 6. - CODSTRAN predicted delect growth. 
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Figure 6. - Continued. 
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Figure i. - Concluded. 
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